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An upper limit to the central density of dark matter haloes 
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ABSTRACT 

We study the growth rates of massive black holes in the centres of galaxies from 
accretion of dark matter from their surrounding haloes. By considering only the ac- 
cretion due to dark matter particles on orbits unbound to the central black hole, 
we obtain a firm lower limit to the resulting accretion rate. We find that a runaway 
accretion regime occurs on a timescale which depends on the three characteristic pa- 
rameters of the problem: the initial mass of the black hole, and the volume density and 
velocity dispersion of the dark matter particles in its vicinity. An analytical treatment 
of the accretion rate yields results implying that for the largest black hole masses 
inferred from QSO studies (> IO^Mq), the runaway regime would be reached on time 
scales which are shorter than the lifetimes of the haloes in question for central dark 
matter densities in excess of 250Mq pc""^. Since reaching runaway accretion would 
strongly distort the host dark matter halo, the inferences of QSO black holes in this 
mass range lead to an upper limit on the central dark matter densities of their host 
haloes of po < 250Mq pc""^. This limit scales inversely with the assumed central black 
hole mass. However, thinking of dark matter profiles as universal across galactic pop- 
ulations, as cosmological studies imply, we obtain a firm upper limit for the central 
density of dark matter in such structures. 
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1 INTRODUCTION 

The question of what is the central density profile of galactic 
dark matter haloes has been much debated in the literature 
over many years. Since the work of Navarro et. al (1997), 
cosmological N-body simulations have consistently agreed in 
yielding central density profiles which can be accurately fit- 
ted by functional forms characterised by centrally divergent 
density cusps. Although the details vary and the innermost 
slope and radius to which said fits are to be trusted are still 
discussed, a broad agreement has been reached in that cos- 
mologically simulated dark haloes exhibit what are termed 
'cuspy density profiles', e.g. Merritt et al. (2006). 

On the other hand, observational inferences of dark halo 
density profiles through rotation curve decomposition, have 
tended to favour those showing density profiles which tend 
to constant values towards the centre. Recent examples in- 
clude de Blok et al. (2008), Kuzio de Naray et al. (2009) and 
Gebhardt & Thomas (2009). The issue is complicated by the 
necessity to estimate the dynamical relevance of the bary- 



onic component, a function of the assumed mass to light 
ratio, the relevance of observational uncertainties such as 
beam smearing, and the importance of non circular motions 
and non centrifugal support of asymmetric drift and hydro- 
dynamic pressure terms (e.g. Valenzuela et al. 2007). 

An interesting independent clue to the puzzle might 
come from the presence of massive black holes in the cen- 
tre of dark matter haloes. The relevance of such black holes 
is well established from their role as the central engines for 
quasars and active galaxies (Rees 1984) as well as in qui- 
escent systems (Kormendy & Richstone 1995). Although 
mostly seen as active nuclei in the distant universe, it is 
commonly believed that all large galaxies host such objects 
in their centres. Recent empirical determinations of central 
QSO black hole masses have established their existence with 
masses in the range 10^ - 10^° Mq at redshifts beyond z~3 
(Kelly et al. 2008, Graham 2008). 

We therefore must conclude that the central regions of 
large dark haloes have coexisted with massive black holes 
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over most of the history of the universe. Given the exis- 
tence of event horizons associated to black holes, and the 
assumption of standard cold dark matter subject only to 
gravitational interactions, it follows that central black holes 
have grown over the history of galactic dark haloes, through 
the capture of dark matter particles. 

The problem of the growth of single central galactic 
black holes through accretion was addressed by Lynden- 
Bell & Rees (1971) and more recently by Gnedin & Primack 
(2004) and Zhao et al. (2002), but only in the accretion of 
particles on capture orbits. While readily acreted, they con- 
stitute only a minor fraction of those available in the dis- 
tribution function of halo particles. Further, once absorbed, 
one has to wait over a comparatively long halo relaxation 
timescale for it to be re-populated with dark matter parti- 
cles. Here we consider only the accretion of unbound parti- 
cles, through the absorption cross section presented by the 
black hole through its event horizon, enhanced by gravita- 
tional focusing. In this case, the accretion is slower at first, 
but proceeds at an ever increasing rate as the growing black 
hole mass leads to an increase in its area, with little accom- 
panying depletion of the overall dark matter halo distribu- 
tion function. This, as accretion does not takes place over 
a highly specific fraction of the distribution function, while 
the black hole mass constitutes only a fraction of the overall 
dark halo mass. 

We find the process to be characterised by the onset of 
a rapid runaway growth phase after a critical timescale. This 
timescale is a function of the mass of the black hole and the 
local density of dark matter. By requiring that the runaway 
phase does not occur, as then the swallowing up of the halo 
by the black hole would seriously distort the former, we can 
obtain upper limits to the maximum allowed density of dark 
matter at the centres of haloes. 



2 CENTRAL BLACK HOLE GROWTH RATES 

In the case of cuspy dark matter profiles, the scales over 
which the density varies by a significant factor, even in the 
central regions, are typically of order 10 pc or above (Merritt 
et al. 2006, Stadel et al. 2009). This is many orders of mag- 
nitude greater than the typical length scales over which the 
accretion onto the central black hole of mass A4 takes place, 
the Schwarzschild radius, iisch = 5 x lQ~''pc(Af/5 x 10^ Mq). 
The range of scales involved, from -Rsch to the scale of the 
dark matter halo, make a full simulation of the whole prob- 
lem unfeasible. To first approximation, we will therefore 
treat the central region of the dark matter halo over which 
the black hole finds itself as one of constant density. 

The effects of central black holes, mostly binary ones, on 
the stellar population of a galaxy have been treated before, 
e.g Quinlan (1996), Merritt & Milosavljevic (2005), Sesana 
et al. (2007). The problem here is different because the total 
mass of the dark halo, and the range of radii covered by the 
dark matter particles are both much larger than equivalent 
quantities in the case of bulge stars affected by central black 
holes. Since the mass of the black hole is still several orders of 
magnitude smaller than the total mass of the dark halo, we 
shall treat the presence of the black hole as a perturbation 
on the distribution function of the dark matter particles. 

In Hernandez & Lee (2008), we calculated the density 



response of a constant density, isothermal dark matter dis- 
tribution, to the presence of a point mass M. A highly lo- 
calised cusp appears, with the density profile changing from 
po to po + pi (r) , where 

^iW = -^Po- (1) 

In the above, a is the velocity dispersion of the halo par- 
ticles, and r is the distance to the mass M. In estimating 
the growth rate of the central black hole, of initial mass Mo, 
we shall begin by assuming that its presence will elicit a re- 
sponse in the otherwise unperturbed dark halo particles, as 
described by eq. (1). In Hernandez & Lee (2008) we showed 
through direct comparison with high resolution N-body sim- 
ulations, that the analytic expression in eq. (1) accurately 
describes the response of a constant density isothermal re- 
gion to the presence of a point mass. A first correction due 
to relativistic effects, the substitution of a Newtonian po- 
tential for the expression of Paczyiiski & Wiita (1980 f], 

will 

result only in the substitution of r — i?sch for the current r 
in the denominator of eq. (1). 

We can now estimate the growth rate of the central 
black hole dimensionally as: 

M = CopAa, (2) 

where p, A and a are a characteristic density, area and 
velocity for the spherical accretion in question, and Co is 
a dimensionless constant which one would expect to be of 
order unity. From the preceding discussion regarding the re- 
action of the halo to the presence of the black hole, and 
from the fact that accretion of unbound particles will occur 
on crossing the event horizon at ii = i?sch, we can estimate 
M from taking p — pi(-Rsch), A = AivRg^^^ and a as the ve- 
locity dispersion of the dark matter particles in the halo. We 
have ignored po in favour of pi in the above considerations, 
as at distances of order Rsch the latter dominates over the 
former by a factor of order {cjcff' . The above yields: 

M = Co8.r^^. (3) 

Gravitational focusing effectively increases the cross 
section of the black hole by trapping particles which would 
otherwise fail to be acreted, and enters into the computation 
of Co. In fact, a fully relativistic calculation for the accretion 
rate of a black hole immersed in an isothermal distribution 
of non-relativistic particles leads to the result: 

M = 16(67r)^''^^!^, (4) 

as derived in Shapiro & Teukolsky (1983), eq. 14.2.26, 
for particles with positive energies. Those on bound orbits 
can be engulfed rapidly in a short initial transient phase, 
and will re-appear as the corresponding phase space is re- 
populated by the distribution function on the very long re- 
laxation timescales of the full halo (Shapiro & Teukolsky 
1983). Here again, by considering only the accretion rate 
of unbound particles, we are confident in having a secure 
lower limit on the accretion rate. Comparing with the di- 
mensional analysis of eq. (3) we see that the result is exact 

^ In the P-W expression a point mass produces a gravitational 
field ^(r) = -~GM/{r - Rsch) instead of the usual $ = -GM/r. 
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for Co = 2(6/7r)i/2_ a factor of less than 3. In what follows 
we shall use the exact result of eq. (4), with cq. (3) serving 
only in allowing a physical interpretation of the relativistic 
result. 

Introducing dimcnsionlcss quantities — a/c, t = t/ta 
and M = M/Mj, where tg = l/(Gpo)^''^ is the free fall 
timescale of the unperturbed background density, and Mj = 
/{G^^^ Pq^) is the Jeans mass of the unperturbed halo, we 
obtain the dimensionless growth rate: 



dM 
dr 



(5) 



We note that in cosmological N-body simulations, the 
distribution function of dark matter particles exhibits a large 
degree of orbital anisotropy and is dominated by highly ra- 
dial orbits (e.g. Ascasibar & Gottlobcr 2008). The increased 
fraction of the dark halo hence available for interaction with 
the central black hole, and the reduced angular momentum 
of the dark matter particles, compared to the isothermal 
case, will all tend to yield faster growth rates than those 
calculated here. The upper limits on central density derived 
below are hence safe upper estimates. More detailed calcula- 
tions accounting for an intrinsically cusped dark halo profile 
with radially dominated distribution functions would yield 
even lower limit densities. 



3 CENTRAL DARK MATTER DENSITY 
LIMITS 

From eq. (4) we obtain for the black hole mass: 

M{t) = 



Moc^a 



c2(T - 16(67r)V2G2poMot' 



(6) 



There is a strong divergence for t tdiv = 

cV/(16(67r)^/^G^poM"o). The time for it to appear decreases 
as the initial black hole mass rises, as the central dark mat- 
ter density increases, and increases as the velocity dispersion 
of the halo particles rises. The divergence is so abrupt that 
the time it takes for the black hole mass to increase by one 
order of magnitude, Tio, is only 9/10 tdiv. The evolution we 
calculate will still be accurate up to f = Tio, as the total 
mass of the dark halo (several times 10^^ Mq for large galac- 
tic haloes), will still be over an order of magnitude larger 
than that of the central black hole, even for the largest ini- 
tial black hole masses considered, of a few times 10^ Mq. We 
shall therefore define: 



Tin = ( ^ 



(7) 



,10/ 16(67r)i/2G2poMo 

as a characteristic timescale after which the accretion pro- 
cess results in substantial dynamical alterations to the over 
all dark halo. With the same dimensionless quantities as de- 
fined previously, we obtain the corresponding expressions: 

Mo 



M{t) 



(8) 



where Ci = 16(67r)i/^ We can now calculate the evolution of 
eq. (6) for any value of the central black hole mass. We begin 
with parameters as appropriate for the largest inferred QSO 
central black holes. Mo = 5 x 10^ M© (Kelly et al. 2008, Gra- 
ham 2008), as this case will lead to the most restrictive dark 
matter density limits. Although the dark haloes of QSOs 



cannot be observationally inferred, given the scalings ob- 
served at low redshift between black hole masses and galac- 
tic properties (e.g. Kormendy & Richstone 1995, Gebhardt 
et al. 2000, Ferrarese & Merritt 2000, Tremaine et al. 2002, 
Gultckin et al. 2009), between black hole masses and halo 
properties (e.g. Volonteri et al. 2003, Bolton et al. 2008, Cro- 
ton 2009, Bandara et al. 2009) and between galactic proper- 
ties and dark halo masses, such as the Tully-Fisher relation 
(Tully & Fisher 1977), it is reasonable to assume that QSOs 
hosting the largest inferred black hole masses will be hosted 
by large galactic haloes. Consequently, we take a large value 
of fT = 200 km/s. Given the usual scaling between cr and the 
flat rotation curve velocity of a galactic halo of Vrot = 2^''^(t, 
this choice corresponds to Vrot = 280 km/s, a value in the 
extreme range for any type of galactic system. Notice that 
as tdiv scales with a, taking a large value for this parame- 
ter will again result in conservative upper limits on the final 
inferred limit central halo densities (see below). Also, given 
the 'inside out' and 'downsizing' aspects of current cosmo- 
logical structure formation models, (e.g. Naab et al. 2009) 
the dynamical stability of the central regions of the most 
meissive galactic systems over the lifetimes we have assumed 
appears reasonable. 

Figure 1 shows the growth of a central black hole as a 
function of time, for our fiducial case with Mo = 5 x 10® M© 
and a = 200 km/s, and a range of values for the assumed 
central dark matter density p()=400, 350, 300, 250, 200, 150 
and 100 AIq pc"'' (from top to bottom, respectively). The 
case where Tio = 10 Gyr corresponds to the middle curve, 
where the central dark matter density is po = 25OM0 pc~^. 
The mass of the black hole increases by a factor of 10 in 
10 Gyr, a conservative estimate of the lifetime of the systems 
in question, namely, dark haloes of QSOs observed at high 
redshift. We see that for central dark matter densities above 
this threshold of 25OM0 pc~^, the mass of the central black 
hole enters the runaway accretion regime and diverges on 
timescalcs shorter than the lifetimes of the systems being 
treated, as given by the three upper curves in Figure 1. On 
the other hand, for values below this threshold, the growth 
of the central black hole is of only a factor of order unity 
over 10 Gyr, as shown by the three lower curves. 

For this particular initial mass. Mo = 5 X 10® Mq, we 
can hence identify 25OM0 pc~^ as a maximum central halo 
dark matter density above which the inferences of black hole 
masses in high redshift QSOs would imply growth rates for 
the central black holes resulting in substantial dynamical 
distortions, leading today not to quiescent black holes in 
the centres of normal galaxies, but to exotic objects dynam- 
ically dominated by extreme super massive black holes. Con- 
sistency arguments of this type can be found e.g. in Gnedin 
& Ostriker (2001), who calibrate the physical parameters of 
self-interacting dark matter by requiring that galactic dark 
haloes should not have evaporated by now into galaxy clus- 
ter dark matter haloes. If we were to take larger values for 
the black hole mass, such as those given by Graham (2008), 
reaching 10^" Mq, or upwards of lO^^M© reported for some 
objects by Kelly et al. (2008), the threshold density we iden- 
tify would go down by a factor of a few. 

In Figure 2 we show the constraints on the central dark 
halo densities, pm, as a function of the assumed central black 
hole mass and for a fixed value of the velocity dispersion, 
cr = 200 km/s. The curves correspond to various values of 
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Figure 1. Evolution of the mass of a central black hole with 
initial mass Mq = 5 x 10^ Mq in a dark halo with dark matter 
particles of isotropic velocity dispersion 200 km/s, for varying 
central region dark matter density: po=400, 350, 300, 250, 200, 
150 and 100 Mq pc~^, top to bottom, respectively. 

Tio. By taking a conservative measure of the lifetimes of 
the systems in question as 10 Gyr, the region containing the 
dotted curves above the thick black line at Tio = 10 Gyr 
is excluded from consistency arguments, while the allowed 
region of parameter space lies in the half plane below it. 
The choice of values for the lifetimes of high redshift QSOs 
larger than the 10 Gyr previously assumed, shifts the max- 
imum central dark matter density values downwards onto 
the various thin solid curves. The most stringent limits ap- 
ply to the highest black hole masses at M ^ 5 x 10^ Mq. 
However, the expectation of universality for the cosmolog- 
ical dark matter density profiles leads one to expect these 
limits will apply to all dark matter haloes. 

We note that once the mass of the central black hole 
grows substantially, processes not included here would begin 
to become relevant, and would invalidate the simple physi- 
cal hypothesis leading to eq. (4). Some include: the adiabatic 
contraction of the dark halo in response to the concentra- 
tion of mass into the central black hole, resulting in higher 
central dark matter densities and hence even higher accre- 
tion rates; the accretion of a fraction of the baryons into the 
central black hole, which is known to occur; or the enhanced 
gravitational focusing of matter of all types into the black 
hole, once the approximation of the black hole mass being 
small compared to the total halo mass which we are working 
under begins to break down. All of these make it reasonable 
to assume that the first corrections to eq. (4) will lead to 
even larger accretion rates, hence leaving our conclusions, 
in terms of limit densities, unchanged. 

Regarding the accretion of baryons and dark halo par- 
ticles, it has been proposed that this can be partly responsi- 
ble for the appearance of the observed scaling laws between 
central black hole masses and bulge and galactic properties, 
e.g., in the analytical work of Zhao et al. (2002) and Gnedin 
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Figure 2. Limit central densities as a function of central black 
hole masses, for a number of values of Tio going from 1 to 15 Gyr, 
every Gyr, from top to bottom respectively. The region above the 
thick black line at Tio = 10 Gyr, with dotted lines, is excluded 
from consistency arguments. 

& Primack (2004), the large scale simulations of Di Matteo 
et al. (2008), or the luminosity function consistency argu- 
ments in Yu & Tremaine (2002) and Hopkins et al. (2007). 
Indeed, Hennawi & Ostriker (2002) constrain the velocity 
dependence of the interaction cross-section of hypothetical 
self-interacting dark matter, by requiring that the accretion 
of dark matter onto central black holes leads to the observed 
M — a relation. Also for the case of self-interacting dark mat- 
ter, Balberg & Shapiro (2002) explore the formation of su- 
permassive black holes through gravothermal core-collapse 
of the central regions of galactic dark haloes. 

Comparing the upper limiting central dark matter den- 
sity of 25OM0 pc~^ with the dynamically inferred structure 
of galactic dark haloes, it is reassuring that when a constant 
density core is used to model observations, the inferred cen- 
tral dark matter densities always lie below this limit, typi- 
cally at ~ IMq pc"'', or below. Recent examples are given 
by Gilmore et al. (2007) for local dwarf spheroidal galajcies, 
and de Blok et al. (2008) for late type galaxies. Hence no 
conflict appears, in that the runaway accretion regime for 
the central black hole will not be reached in 10 Gyr for any 
directly inferred values of the central dark matter density, 
for any inferred central black hole masses. 

From the point of view of cuspy dark matter haloes, 
the limits we derive here establish an inner boundary, ex- 
terior to which the globally fitted centrally divergent dark 
halo profiles can be valid. At smaller radii, this solution 
must be modified to avoid the divergent black hole growth 
rates found here. Although our high limit densities are not 
reached by current cosmological N-body simulations, which 
typically stop at volume densities of order IMq pc"'' at 
their resolution limit, even for the most recent highest reso- 
lution experiments (Cuesta et al. 2008, Stadel et al. 2009), 
the logarithmic slopes in these regions are still such that 
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volume densities above the limits we derive here would be 
reached orders of magnitude in the radial coordinate before 
reaching the scale of the super massive central black holes, 
even for central black hole masses below the upper ranges of 
these values. Studies of the origin of cuspy cosmological den- 
sity profiles within the secondary infall scenario have traced 
the cusp to the close to scale free initial perturbation spec- 
trum (e.g. Williams et al. 2004, Salvador-Sole et al. 2005, Del 
Popolo 2009). Such studies explain the steep negative loga- 
rithmic slopes of density profiles from N-body simulations, 
predicting them to extend into the very centres. Although 
these considerations apply only to cosmological dark matter 
haloes in the absence of any central black holes, the inclu- 
sion of single central black holes, currently beyond the reach 
of fully self consistent simulations, will result in even steeper 
central dark matter profiles, strengthening the consistency 
arguments made here. A solution might include dark mat- 
ter physics not ordinarily considered, such as self-interacting 
dark matter (e.g. Firmani et al. 2000, de la Macorra 2009), 
warm dark matter, or changes to the initial fiuctuation spec- 
trum (e.g. Alam et al. 2002). 



4 CONCLUSIONS 

We study the mass growth rates of central black holes 
through accretion of dark matter particles on orbits un- 
bound to the central black hole. As the black hole mass 
grows, a runaway accretion regime ensues. Requiring that 
no such runaway regime has been reached over lifetimes of 
galactic dark haloes of 10 Gyr leads to the identification of 
critical upper limits for the density of dark matter. These 
limits scale proportionally to the assumed value of the dark 
matter velocity dispersion, and inversely proportionally to 
the assumed value of the central black hole mass. 

For the largest black hole masses inferred by QSO stud- 
ies of 5 X 10^ Mq, central region dark matter densities larger 
than po = 25OM0 pc~'^ are excluded. These limits suggest 
dark halo density structures are characterised by constant 
density central regions, rather than divergent cuspy profiles. 
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